AD-A246  963 


E 


la  REPORT  SECURITY  CLASSIFICATION 


a.  SECURITY  CLASSIFICATION  AUT 


2d  DECLASSiFICATiON  i  DOWNGRA 


TATI  ON  PAGE 


lb  RESTRICTIVE  MARKINGS 


3  Distribution /availability  of  report 


Approved  for  public  release; 
distribution  is  unlimited. 


Form  Approved 
OMB  No  0704-01 86 


4  PERFORMING  ORGANIZATION  REPORT  NUMBER(S 


6a.  NAME  OF  PERFORMING  ORGANIZATION 

Georgia  Univ  Research 
Foundation  Inc 

6b  OFFICE  SYMBOL 
(If  applicable) 

6c.  ADORESS  (City,  State,  and  ZIP  Code) 

Athens,  GA  30602 

8a.  NAME  OF  FUNDING/SPONSORING 

8b  OFFICE  SYMBOL 

ORGANIZATION 

(If  applicable) 

AFOSR 

NC 

8c.  ADDRESS  (City,  State,  and  ZIP  Code) 

Building  410,  Bolling  AFB  I 

)C 

20332-6448 

AFOSR/NC 


7b  ADDRESS  (City.  Stats.  and  ZIP  Code) 

Building  410,  Bolling  AFB  DC 
20332-6448 


AFOSR-88-0167 


10  SOURCE  OF  FUNDING  NUMBERS 


PROJECT 

TASK 

NO 

NO 

2303 

B3 

WORK  UNIT 
ACCESSION  NO 


11.  TITLE  (Include  Security  Classification ) 

(U)  FUNDAMENTAL  STUDIES  OF  CARBON,  NH,  AND  OXYGEN  RINGS  AND  OTHER  HIGH 
ENERGY  DENSITY  MOLECULAR  SYSTEMS 


12  PERSONAL  AUTHOR(S) 


13a.  TYPE  OF  REPORT 

FINAL 


16.  SUPPLEMENTARY  NOTATION 


13b  TIME  COVERED  14.  DATE  OF  REPORT  (Year,  Month,  Day)  15.  PAGE  COUNT 

from  to  in-31 -91  December  30,  1991 


COSATI  CODES 


FIELD  GROUP 


SU8-GROUP 


18  SUBJECT  TERMS  ( Continue  on  reverse  if  necessary  and  identify  by  block  number) 

Ab  Initio,  computational  chemistry,  quantum  chemistry, 
theoretical  chemistry,  propellants 


19  A8STRACT  ( Continue  on  reverse  if  necessary  and  identify  by  block  number) 

The  object  of  this  research  is  to  characterize  the  molecular  structures, 
energetics,  spectroscopic  properties,  and  elementary  chemical  reactions  of 
the  oxygen  ring  molecules  O4  through  Oj£  and  related  species  including  (NH)n 
and  Cn.  The  approach  used  will  exploit  recent  developments  in  ab  initio  molecular 
quantum  mechanics. 


20  DISTRIBUTION /AVAILABILITY  OF  ABSTRACT  21  A8STRACT  SECURITY  CLASSIFICATION 

□  unclassified/unlimited  □  same  as  RPT  Q  otic  USERS  UNCLASSIFIED 


22a  NAME  OF  RESPONSIBLE  individual  22b  'E-EP^ONE  (Include  Area  Code)  22 c  OFFICE  SYMBOL 

Dr.  Mark  S.  Gordon  (202)  767-4963  AFOSR/NC 


DO  Form  1473,  JUN  86 


Previous  editions  are  obsolete 


SECURITY  CLASSIFICATION  QF  *hi$  ;aQE 

UNCLASSIFIED 


r 


t 


FINAL  TECHNICAL  REPORT 


Ac«««alda  F*r 

NTI'J  QjUAl 
»T4C  tib  □ 

1  Juatlf  leatloo. _ 


Dl3trtlMttl*u/ 


AvaUatUlty  Co««e 

£**•!!  a»4/«r 


Air  Force  Office  of  Scientific  Research 
Grant  AFOSR-88-0167 


Diet 

|  Spe«lal 

(V\ 

’’Fundamental  Studies  of  Carbon,  NH,  and  Oxygen  Rings 
and  Other  High  Energy  Density  Molecular  Systems” 


I.  Summary 

The  development  of  efficient  and  safe  conventional  (i.e.,  non-nuclear)  propellants 
and/or  fuels  is  a  goal  of  obvious  technological  significance.  A  desirable  quality  of  such  a 
propellant  is  clearly  a  high  ratio  of  energy  release  to  mass.  The  present  research  rests  on 
a  simple,  but  previously  (i.e.,  before  our  work  beginning  in  1987)  unrecognized,  analogy 
between  oxygen  and  sulfur.  Our  AFOSR  supported  research  of  the  past  four  years  has 
shown  that  the  proposed  oxygen  ring  systems  are  sufficiently  promising  to  warrant 
further  consideration.  In  fact,  several  experimental  groups  are  now  attempting  to  prepare 
and  characterize  oxygen  rings  in  the  laboratory.  Thus  a  major  emphasis  in  future  phases 
of  the  HEDM  (high  energy  density  molecules)  theoretical  research  is  to  work  closely 
with  Air  Force  supported  experimental  efforts  to  this  end.  Our  work  has  also  encouraged 
other  theoretical  groups  to  extend  our  S  — >  O  analogy  to  P  N,  with  the  recent  work  of 
Lee  and  Rice  on  tetrahedral  N4  being  a  particularly  beautiful  example. 

The  potential  attractiveness  of  oxygen  rings  follows  from  a  number  of  considera- 
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(a)  Our  theoretical  studies  to  date  show  that  the  cyclic  On  systems  are  definitely  high 
energy  density  materials.  This  may  be  seen  in  Figure  1,  kindly  provided  to  us  by 
Dr.  Harvey  Michels  based  on  our  predicted  energetics  for  Og; 

(b)  As  cryogenic  solid  propellants,  oxygen  rings  may  plausibly  be  expected  to  form  a 
metastable  condensed  phase  by  analogy  with  the  valence  isolectronic  sulfur  rings, 
which  of  course  are  very  stable  solids  under  normal  conditions; 

(c)  The  way  to  use  oxygen  rings  as  a  rocket  fuel  is  obviously  analogous  to  that  for  nor¬ 
mal  molecular  oxygen,  i.e.,  O2.  Furthermore,  burning  of  On  with  molecular  hydro¬ 
gen  yields  water  as  the  only  new  molecular  product.  The  absence  of  hazardous 
products  (as  opposed  to  other  possible  HEDM  materials)  of  combustion  is  a  strong 
argument  in  favor  of  oxygen  rings. 

Let  us  summarize  the  qualitative  arguments  concerning  the  energy  content  of  the 
metastable  On  rings.  First  one  notes  that  the  valence-isoelectronic  cyclic  S8  contains  no 
energy  in  this  sense.  Specifically,  gaseous  Sg  lies  12.2  kcal  per  mole  of  sulfur  atoms 
below  four  gas  phase  diatomic  S2  molecules.  Can  one  guess  this  result  using  the  dissoci¬ 
ation  energy  of  S2  (D0  =  4.37  eV  =  100.8  kcal/mole)  and  some  reasonably  standard  S-S 
single-bond  dissociation  energy,  say  54  kcal/mole?  The  answer  to  this  question  is  a 
qualified  "yes".  Using  this  simple  model,  one  predicts  Sg  to  lie  (54-50.4)  =  3.6 
kcal/mole  below  four  S2  molecules  on  a  per  atom  basis.  By  increasing  the  S-S  bond 
energy  from  the  standard  54  kcal  to  63  kcal,  the  known  experimental  energy  difference 
(+97.7  kcal)  for 


Sg(g)  4  S2(g) 


(1) 


H2/08  Vacuum  Specific  Impulse 

Chamber  Pressure  =  1000  psia 


Oxidizer/Fuel  Weight  Ratio 

Figure  1 .  Specific  impulse  predictions  for  cyclic  08. 
Details  are  given  in  the  text. 
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is  precisely  reproduced.  This  adjustment  reflects  the  fact  that  the  S-S  bonds  in  Sg  are 
stronger  than  those  in  organosulfiir  compounds  such  as  CH3SSCH3. 

The  comparable  oxygen  thermodynamic  data  make  it  immediately  obvious  why 
oxygen  rings  should  be  high  energy  density  materials.  First  the  standard  0-0  bond  dis¬ 
sociation  strength  is  35  kcal/mole,  much  weaker  than  the  54  kcal  for  S-S  bonds. 
Secondly,  the  dissociation  energy  of  diatomic  O2  is  D0  =  5.12  eV  =  118.0  kcal/mole, 
much  stronger  than  the  100.8  kcal  for  S2.  Thus  the  estimate  for  the  energy  stored  in  the 
"generic"  oxygen  ring  is  (35-59)  =  -24  kcal  per  mole  of  oxygen  atoms.  That  is,  the  dis¬ 
sociation  process 

Oa  (n  even)  ->  y  02  (g)  (2) 

is  estimated  to  be  exothermic  by  24  kcal  per  mole  of  oxygen  atoms.  Alternately,  48  kcal 
of  energy  is  released  per  mole  of  02  molecules.  A  larger  energy  release  is  expected  for 
the  smaller  rings,  specifically  04  and  0&,  which  presumably  have  smaller  average  0-0 
bond  energies  due  to  ring  strain. 

A  principal  accomplishment  of  the  current  AFOSR  grant  has  been  the  prediction  of 
the  molecular  structures  and  energetics  of  04,  Og,  and  O12.  References  to  our  published 
work  on  these  systems  are  given  below.  A  careful,  high  level  theoretical  study  of  06  is 
currently  underway.  Figures  2  and  3  present  our  theoretical  structures  for  Og  and  0|2- 
We  are  confident  that  these  structures  will  be  reliable  should  the  oxygen  rings  prove  sus¬ 
ceptible  to  detailed  molecular  structure  investigations  in  the  laboratory. 

What  have  we  learned  about  the  energetics  of  the  oxygen  rings?  Let  us  first  focus 
on  the  results  for  O12.  The  valence  isoelectronic  sulfur  molecule  S12  at  the  DZ  SCF 
level  of  theory  is  predicted  to  lie  21.9  kcal/mole  below  six  S2  molecules.  This  result  may 
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readily  be  translated  into  the  prediction  that  S12  lies  21.9/12  =  1.8  kcal/mole  below 
separated  S2  molecules  on  a  per  atom  basis.  The  analogous  experimental  value  for  gase¬ 
ous  Sg  is  12.2  kcal/mole  per  S  atom.  It  is  of  course  anticipated  that  DZ  SCF  theory  will 
do  a  better  job  for  6  S2  than  for  Si2,  and  an  error  of  -10  kcal/mole/atom  seems  to  us  to 
be  perhaps  less  than  might  be  expected.  At  the  DZ  SCF  level  012  lies  20.9  kcal/mole 
above  six  02  molecules  on  a  per  atom  basis.  With  the  larger  basis  set  DZ+P  SCF  method 
Ol2  lies  21.6  kcal/mole/atom  above  six  separated  diatomic  oxygen  molecules.  Thus  the 
addition  of  polarization  functions  (d  functions  on  each  oxygen  atom)  lowers  the  energy 
of  6  02  molecules  somewhat  more  than  that  of  0I2. 

Analogous  to  S|2,  one  certainly  expects  that  higher  levels  of  theory  (larger  basis 
sets,  and  especially  explicit  treatment  of  electron  correlation)  will  lower  the  energy  of 
012  relative  to  6  02.  Second-order  perturbation  theory  supports  this  view.  With  the  DZ 
basis  set,  the  MP2  energy  difference  between  012  and  6  02  is  16.1  kcal/mole  on  a  per  O 
atom  standard.  These  results  suggest  that  0J2  may  be  significantly  more  stable  than  our 
earlier  back-of-the-envelope  calculation  (based  on  standard  bond  energies),  which  as  dis¬ 
cussed  above  indicated  that  generic  cyclic  0D  molecules  might  store  24  kcal  per  mole  of 
oxygen  atoms  relative  to  the  separated  02  molecules. 

Turning  to  the  eight-membered  rings,  reported  in  Table  I  are  the  energies  of  Og  and 
Sg  relative  to  the  dissociation  limits  4  02  and  4  S2.  Since  experimental  data  is  available 
for  S<  these  results  will  be  discussed  first.  At  the  DZ  SCF  and  DZ+P  SCF  levels  of 
theory,  Sg  lies  13.0  and  33.9  kcal/mole,  respectively,  below  four  separated  S2  molecules. 
The  latter  (DZ+P  SCF)  result  may  readily  be  translated  into  the  prediction  that  Sg  lies 
33.9/8  =4.2  kcal/mole  below  separated  S2  molecules  on  a  per  atom  basis.  The  analo¬ 
gous  experimental  value  for  gaseous  S8  is  12.2  kcal/mole  perS  atom. 
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Assuming  the  DZ+P  SCF  equilibrium  geometry  for  S8  (Figure  3),  its  dissociation 
energy  has  been  predicted  using  second-order  perturbation  theory  (MP2).  However,  for 
an  open-shell  system  like  the  S2  molecule  3  Eg  ground  state,  the  use  of  unrestricted 
Mriller-Plesset  second-order  perturbation  theory  is  most  straightforward.  At  the  S2  RHF 
equilibrium  geometry  (DZ+P  basis  set,  re  =  1.881  A)  the  UMP2  total  energy  is 
-795.25545  hartrees.  The  comparable  MP2  energy  for  closed  shell  Sg  is  -3181.11839 
hartrees.  Thus  Sg  is  predicted  to  lie  60.6  kcal/mole  below  4  S2  at  this  level  of  theory. 
On  a  per  atom  basis  Sg  lies  7.6  kcal/mole  below  4  S2,  in  reasonable  agreement  with  the 
experimental  value  12.2  kcal/mole. 

At  the  DZ  SCF  level  Og  lies  21.3  kcal/mole  above  four  O2  molecules  on  a  per  atom 
basis.  With  the  larger  DZ+P  basis  set  Og  lies  21.6  kcal/mole/atom  above  four  separated 
diatomic  oxygen  molecules  within  the  SCF  approximation.  Thus  the  addition  of  polari¬ 
zation  functions  (d  functions  on  each  oxygen  atom)  lowers  the  energy  of  4  02  slightly 
more  than  that  of  Og. 

Assuming  DZ+P  SCF  geometries,  the  DZ+P  MP2  energies  for  02  and  Og  are 
-149.97388  and  -599.68739  hartrees,  respectively.  Thus  Og  lies  130.6  kcal/mole  above  4 
O2  or  16.3  kcal/mole  higher  on  a  per  atom  basis.  Hence  the  effect  of  electron  correlation 
is  to  significantly  lower  Og  relative  to  the  separated  oxygen  molecules. 

More  reliable  yet  should  be  the  comparison  of  the  energy  of  DZ+P  MP2  geometry 
optimized  Og  with  that  of  four  comparable  O2  molecules.  The  O2  geometry  optimiza¬ 
tion  is  carried  out  at  the  DZ+P  UMP2  level,  yielding  a  bond  distance  re  =  1.253  A  and 
total  energy  E  =  -149.97905  hartrees.  Combined  with  the  Og  total  energy,  one  predicts 
that  cyclooctaoxygen  lies  123.5  kcal/mole  above  four  infinitely  separated  O2  molecules. 
Thus  geometry  optimization  is  energetically  much  more  important  for  Og  with  correlated 
methods  than  is  the  case  for  O2.  On  a  per  atom  basis  the  DZ+P  MP2 
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dissociation  energy  for  Og  is  15.4  kcal/mole. 

A  number  of  significant  relationships  between  the  Og  energetic  predictions  and 
those  for  0,2  may  be  noted.  First,  the  higher  level  theoretical  methods  used  are  broadly 
in  agreement  with  the  simpler  methods  used  for  O12.  For  0,2  the  only  correlated  method 
used  was  DZ  MP2  (using  DZ  SCF  geometries)  and  those  results  were  presented  with 
great  caution.  However,  the  analogous  DZ  MP2  and  DZ+P  MP2  dissociation  energies 
for  Og  are  15.7  and  16.3  kcal/mole/atom,  respectively.  The  close  agreement  between  the 
two  methods  suggests  that  our  final  0,2  energetic  predictions  may  be  much  more  reliable 
than  could  reasonably  have  been  anticipated. 

Secondly  the  dissociation  energies  of  Og  and  O12  appear  to  be  very  similar  on  a  per 
atom  basis.  This  may  be  seen  from  the  following  array  of  dissociation  energies: 


Og 

012 

DZ  SCF 

21.3 

20.9 

DZ+P  SCF 

21.6 

21.6 

DZ  MP2 

15.7 

16.1 

Both  Og  and  012  should  be  relatively  free  of  bond  angle  and  dihedral  angle  strain  and 
their  comparable  energetics  mirror  those  presumed  for  Sg  and  S,2  based  on  the  latter’s 
stability  (the  precise  thermochemistry  of  gaseous  S,2  does  not  appear  to  be  established). 
Thus  one  expects  Og  and  012  to  be  perhaps  the  most  readily  synthesizable  of  the  oxygen 
rings. 

Among  the  oxygen  rings,  the  energetics  have  been  most  reliably  established  for 
what  is  apparently  the  least  stable  system,  namely  cyclotetraoxygen  (Figure  4).  Thus  we 
are  now  able  to  make  a  reliable  prediction  of  the  heat  of  formation  of  cyclic  O4. 

Using  DZ+P  CCSD  energies,  we  can  estimate  the  enthalpy  change  for  converting 
four  moles  of  O3  to  three  moles  of  04.  The  electronic  contribution  to  this  enthalpy 


cyclic 
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change  is  112.6  kcal/3  moles.  The  rotational  and  vibrational  contributions  to  this 
enthalpy  change  total  -1.5  kcal/3  moles.  The  experimental  enthalpy  of  formation  for  03 
is  34.1  kcal/mole  or  136.4  kcal/4  moles.  The  enthalpy  of  formation  of  04  can  be 
estimated  by  adding  the  enthalpy  of  formation  for  four  moles  of  03  to  the  enthalpy 
change  for  converting  four  moles  O?  to  three  moles  04.  The  enthalpy  of  formation  of  04 
is  therefore  approximately  247.5  kcal/3  moles  or  82.5  kcal/mole.  The  latter  prediction  is 
17.5  kcal/mole  lower  than  our  previous  estimate  based  on  the  DZ+P  C1SD  level  of 
theory.  Some  of  the  difference  comes  from  the  use  of  the  experimental  heat  of  forma¬ 
tion  for  ozone.  Lf  this  had  been  done  for  the  CISD  estimate,  the  result  would  have  been 
94.3  kcal/mol.  Because  of  the  size  consistency  of  the  CCSD  method,  the  CCSD  estimate 
is  deemed  to  be  more  accurate  than  CISD.  Thus  the  energetics  of  the  04  system  at  equili¬ 
brium  seem  to  be  quite  reliably  established. 
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